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Abstract 

In this paper, we propose a simple rule that generates scale-free small-world net- 
works with tunable assortative coefficient. These networks are constructed by two- 
stage adding process for each new node. The model can reproduce scale-free degree 
distributions and small- world effect. The simulation results are consistent with the 
theoretical predictions approximately. Interestingly, we obtain the nontrivial clus- 
tering coefficient C and tunable degree assortativity r by adjusting the parameter: 
the preferential exponent (3. The model can unify the characterization of both as- 
sortative and disassortative networks. 
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1 Introduction 



In the past few years, no issues in the area of network researching attract 
more scientists than the ones related to the real networks, such as the Internet, 
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the World-Wide Web, the social networks, the scientific collaboration and so 
on [1,2,3,4,5,6]. Recent works on the complex networks have been driven by the 
empirical properties of real-world networks and the studies on network dynam- 
ics [7,8,9,10,11,12,13,14,15,16,17,18,19,20,21,22,23,24,25,26,27,28,29]. Many em- 
pirical evidences indicate that the networks in various fields have some com- 
mon topology characteristics. They have a small average distance like random 
graphs, a large clustering coefficient and power-law degree distribution [1,2], 
which are called the small-world and scale-free characteristics. The other char- 
acteristic is that the social networks are assortative while almost all biological 
and technological networks are opposite. The networks with high clustering 
and small average distance are the small-world model of Watts and Strogatz 
(WS)[1], while the networks with power- law degree distribution are the scale- 
free network model of Barabdsi and Albert (BA) [2] . The BA model is a pio- 
neering work in the studies on networks, which suggests that the growth and 
preferential attachment are two main self-organization mechanisms. Although 
BA model can generate the power-law degree distributions, its assortative co- 
efficient r equals to zero in the limit of large size thus fail to reproduce the 
disassortative property that extensively exists in the real-world networks. Re- 
cently, some models that can generate either assortative or disassortative net- 
works have been reported [30,31,32,33,34,35]. Wang et al. presented a mutual 
attraction model for both assortative and disassortative weighted networks. 
The model found that the initial attraction A of the newly added nodes may 
contribute to the difference of the assortative and disassortative networks [34]. 
Liu et al. [36] proposed a self-learning mutual selection model for weighted 
networks, which demonstrated that the self-learning probability p may be the 
reason why the social networks are assortative and the technological networks 
are disassortative. However, one should not expect the existence of a omnipo- 
tent model that can completely illuminate the underlying mechanisms for the 
emergence of disassortative property in various network systems. In this pa- 
per, beside the previous studies, we exhibit an alternative model that can 
generate scale- free small- world networks with tunable assortative coefficient, 
which may shed some light in finding the possible explanations to the different 
evolution mechanisms between assortative and disassortative networks. 

Dorogovtsev et. al [37] proposed a simple model of scale-free growing net- 
works . In this model, a new node is added to the network at each time step, 
which connects to both ends of a randomly chosen link undirected. The model 
can be equally described by the process that the newly added node connect 
to node i preferentially, then select a neighbor node of the node i randomly. 
Holme and Kim [38] proposed a model to generate growing scale-free networks 
with tunable clustering. The model introduced an additional step to get the 
trial information and demonstrated that the average number of trial informa- 
tion controls the clustering coefficient of the network. It should be noticed 
that the newly added node connect to the first node i preferentially, while 
connect to the neighbor node of the first node % randomly. In this paper, we 
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will propose a growing scale-free network model with tunable assortative co- 
efficient. Inspired by the above two models, the new node is added into the 
network by two steps. In the first step, the newly added node connects to the 
existing nodes i preferentially. In the second step, this node selects a neighbor 
node s of the node i with probability fcf / X^er, kf, where (3 is the parameter 
named preferential exponent and Tj is the neighbor node set of node %. This 
model will be equal to the Holme-Kim model[37] when (3 = 0, and the MRGN 
model[35] when j3 — 1 Specifically, the model can generate a nontrivial clus- 
tering property and tunable assortativity coefficient. Therefore, one may find 
explanations to various real-world networks by our microscopic mechanisms. 



2 Construction of the Model 

Our model is defined as following. 

(1) Initial condition: The model starts with m connected nodes. 

(2) Growth: At each time step, one new node v with m edges is added at 
every time step. Time t is identified as the number of time steps. 

(3) The first step: Each edge of v is then attached to an existing node with 
the probability proportional to its degree, i.e., the probability for a node i 
to be attached to v is 
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(4) The second step: If an edge between v and w was added in the first step, 
then add one more edge from v to a randomly chosen neighbor s of w with 
probability P t according to the following probability 



k 13 

P, = -^J, (2) 



If there remains no pair to connect, i.e., if all neighbors of w were always 
connected to v, do the first step instead. 
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3 Characteristics of the Model 



3.1 Degree distribution 



The degree distribution is one of the most important statistical character- 
istics of networks. Since some real-world networks are scale-free, whether the 
network is of the power-law degree distribution is a criterion to judge the va- 
lidity of the model. By adopting the mean-field theory, the degree evolution 
of individual node can be described as 

flk 

= P(i) + E P(i\j)PV), (3) 



where P(i) denotes the probability that the node i with degree ki is selected 
at the first step, P(i\j) denotes the conditional probability that node % is a 
neighbor of node j with degree kj which has been selected at the first step. 

According to the preferential attachment mechanism of the first step, one 
has 

p(i) = A < 4 » 



The conditional probability P(i\j) can be calculated by 



According to the second step, one has that 
dkj _ h / fcf kj \ 

dt Eik ftr^Es^klEik'- K> 



If (3 — 0, we get that 

dki ki f 1 kj \ 2/uj . . 

+ E (tT^T-J = ~7T- (7) 



dt Ei h f-£. v kj Ei k ' Ei h 



Then we can get that P(k) ~ k 3 , which has been proved by Holme and Kim 
[38]. If /? = 1, the following formula can be obtained under the assumption 
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that the present network is non-assortative. 

dki ki ^ ki k\ 2ki 

dt J2j kj (k)ki kj J2j kj 

We can obtain that the degree distribution p(k) ~ k~ r obeys the power-law 
and the exponent 7 = 3. The numerical results are demonstrated in Fig. 
1. From Fig. 1, we can get that the exponents 7 of the degree distribution 
are around -3 when (3 = —2, —1,0, 1,2. When (3 > 0, the exponent 7 would 
increase slightly as the j3 increases. 



q — 1 1 1 1 1 1 1 1 — 1 — 1 1 1 1 1 1 1 1 — 1 — 1 1 1 1 1 1 1 1 — 1 — 1 1 1 1 1 1 1 1 — 1 — 1 1 1 1 1. 
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k 

Fig. 1. (Color Online) Degree distribution of the present network with N = 100000 
nodes when m = 3 and Pt = 0.3. In this figure, p(k) denotes the probability of 
nodes with degree k in the network. The power-law degree distribution exponents 7 
of the four probability density function are 7/3=2 = 3.11 ± 0.05, 7/3=1 = 3.11 ± 0.05, 
7/3=0 = 2.93 ± 0.04, 7/3=^ = 2.96 ± 0.05 and 7/3=_ 2 = 2.95 ± 0.04. 

3.2 Average distance 

The average distance is also one of the most important parameters to mea- 
sure the efficiency of communication networks, which is defined as the mean 
distance over all pairs of nodes. The average distance plays a significant role 
in measuring the transmission delay. Firstly, we give the following lemma [39]. 
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Fig. 2. (Color Online) The dependence between the average distance D and the order 
N of the present network, when (3 = 1 and m = 4. One can see that L increases 
very slowly as IniV increases. The inset exhibits the curve where L is considered as 
a function of N. All the data are obtained by 10 independent simulations. 

Lemma 1 For any two nodes i and j , each shortest path from i to j does 
not pass through any nodes k satisfying that k > max{i, j}. 

Proof. Denote the shortest path from the node i to j of length n + 1 by 
i — > x 1 — > x 2 ■ ■ ■ — > x n — > j(SPij), where n > 0. Suppose that x k = 
maxjx 1 , x 2 , ■ ■ ■ , x n }, if x k < max{i,j}, then the conclusion is true. If x k > 
max{i,j}, denote the youngest node of SPij by k. Denote the subpath pass- 
ing through node k by I — > k — > r, where the node I and r are the neighbors 
of node k, then we can prove that node I and r are connected. The shortest 
path SPij passes from the node I to r directly, which is conflicted with the 
hypothesis. 

Let d(i,j) represent the distance between node % and j and cr(iV) as the total 
distance, i.e., cr(N) = Ei<i<i<jv j). The average distance of the present 
model with order N, denoted by L(N), is defined as following 

L{N) = M") (9) 

According to Lemma 1, the newly added node will not affect the distance 



6 



between the existing ones. Hence we have 



N 



a(N + l) = a(N) + J2 d (^ N + ( 10 ) 



Assume that the (N + l)th node is added to the edge E yW2 , then Equ. (10) 
can be rewritten as 

a(N + 1) = a(N) + N + f] D(i, y), (11) 

i=i 



where D(i,y) = mm{d(i,yi),d(i,y2)}- Denote y as the edge connected the 
node y± and y 2 continuously, then we have the following equation 

a(N + 1) = a(N) + iV + ]T d(i, y), (12) 

i=A 



where the node set A = {1, 2, • • • , N} — {y 1 ,y 2 } has (N — 2) members. The sum 
J2i=\d(i,y) can be considered as the distance from each node of the network 
to node y in the present model with order N — 2. Approximately, the sum 
J2i=\ d(h y) is equal to L(N — 2). Hence we have 

J2d(i,y)*(N-2)L(N-2). (13) 

i=A 



Because the average distance L(N) increases monotonously with N, this yields 

(N-2)L(N-2) = (N-2), MN f 2) ; < (14) 
V ; V ; V ; (A^-2)(AT-3) JV-1 V ; 



Then we can obtain the inequality 

<7 (JV + 1)< <7 (JV) + JV+?^. (15) 



Enlarge cr(N), then the upper bound of the increasing tendency of cr(N) reads 



This leads to the following solution 

a(N) = (N — l) 2 log(N - 1) + d(N - l) 2 — (N — 1). (17) 
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Fig. 3. (Color Online) The scale of C with various (3 and Pt when m = 3. The data 
are averaged over 10 independent runs of network size of N = 6000. 



The numerical results are demonstrated in Fig. 2. 



3. 3 Clustering property 

The small-world characteristic consists of two properties: large clustering 
coefficient and small average distance. The clustering coefficient, denoted by 
C, is defined as C = YliLx %, where Cj is the local clustering coefficient for 
node i. C{ is 



where E(i) is the number of edges in the neighbor set of the node i, and hi 
is the degree of node i. When the node % is added to the network, it is of 
degree m + mP t and E(i) = mP t . If a new node is added to be a neighbor of i 
at some time step, E(i) will increase by mP t since the newly added node will 
connect with one of the neighbors of the node % with probability Pt. Therefore, 
in terms of hi, the expression of E(i) can be written as 



2E(i) 



(18) 



hi(hi 1) 



E(i) = mP t + P t [ki - (m + mP t )]. 



(19) 



t 1 1 1 — i — i — r-| 1 1 1 1 1 — i — i — r-| r 

10 100 

k 



Fig. 4. (Color Online) C(k) vs k to various f3 and when m = 3 and Pt = 0.3. The 
data are averaged over 10 independent runs of network size of N = 6000. 



Hence, we have that 

( k ■ — m — mP.W / Tti P.2 P, — m P? \ 

(20) 



c = 2[mP t + P t (kj - m - mP t )} = ^/mP^ P t - mP t 
k>i{k"i 1) k% k% 1 



This expression indicates that the local clustering scales as C(k) ~ where 
C(k) denotes the average clustering coefficient value of nodes with degree k. It 
is interesting that a similar scaling has been observed in many artificial models 
[38,39,40,41] and several real- world networks [42]. The degree-dependent aver- 
age clustering coefficient C(k) has been demonstrated in Fig. 4. Consequently, 
we have 




Since the degree distribution is P(k) = ak 3 , where k = k min , • • • , fc max - The 
constant a satisfies the normalization equation 

E ak ~ 3 = !> ( 22 ) 

k ■ 
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one can get that a = 2/c^ in . The average clustering coefficient C can be rewrit- 
ten as 



r _ 2_ V U ( NP(k)mP? NP(k){P t -mP?y 

N ^femin V k I" k-1 y ^ 

~ 2 Efc (ak^mP? + ak- 3 (P t - mP t 2 )/(k - 1)) 



The numerical results are demonstrated in Fig.3. From figure 3, we can get that 
if f3 < 0, the numerical results are consistent with the theoretical predictions 
approximately, while if (3 > 0, the fluctuations emerges. The departure from 
analysis results is observed, which may attribute to the fluctuations of the 
power-law exponent of degree distribution. It is also helpful to compare the 
present method with previous analysis approaches on clustering coefficient for 
Holme-Kim model [43,44]. 



3.4 Assortative coefficient 



The assortative coefficient r can be calculated from 



where ji, ki are the degrees of the vertices of the ith edge, for i — 1, 2, • • • , M 
[45,46]. 

From Fig. 5, we can find that when (3 > the assortative coefficient r 
increases with the probability Pt, while r decreases with the probability Pt 
when P < 0. As/3 = 0,r equals to zero approximately. 



4 Conclusion and Discussion 



In this paper, we propose a simple rule that generates scale-free small-world 
networks with tunable assortative coefficient. The inspiration of this model 
is to introduce the parameter f3 to Holme-Kim model. The simulation results 
are consistent with the theoretical predictions approximately. Interestingly, we 
obtain the nontrivial clustering coefficient C and tunable degree assortativity 
r, depending on the parameters (3. The model can unify the characterization of 
both assortative and disassortative networks. Specially, studying the degree- 
dependent average clustering coefficient C(k) also provides us with a better 
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Fig. 5. (Color Online) The scale of r with various (3 and Pt when m = 1. The data 
are averaged over 10 independent runs of network size of N = 6000. 

description of the hierarchies and organizational architecture of weighted net- 
works. Our model may be conducive to future understanding or characterizing 
real-world networks. 

This work has been supported by the Chinese Natural Science Foundation 
of China under Grant Nos. 70431001, 70271046 and 70471033. 
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